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We have performed AC susceptibility and DC magnetic relaxation measurements on the spin ice system Dy2Ti207 
down to 0.08 K. The relaxation time of the magnetization has been estimated below 2 K down to 0.08 K. The spin 
dynamics of Dy2Ti207 is well described by using two relaxation times (ts (short time) and tl (long time)). Both ts and 
Tl increase on cooling. Assuming the Arrhenius law in the temperature range 0.5-1 K, we obtained an energy barrier 
of 9 K. Below 0.5 K, both ts and tl show a clear deviation from the thermal activated dynamics toward temperature 
independent relaxation, suggesting a quantum dynamics. 
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The pyrochlore oxide Dy 2X1267 is considered to be a 
typical example of a spin ice compound. In this com- 
pound, Dy^^ ions occupy the vertices of corner-shared tetra- 
hedra. The magnetic moment of Dy^^ has a strong single- 
site anisotropy along the local (111) axes, which is caused 
by crystalline electric field eff'ects (CEF) and gives rise to a 
Kramers ground state doublet. The strong dipolar interaction 
between Dy^^ ions associated with the weak antiferromag- 
netic superexchange interaction results in an eff'ective ferro- 
magnetic nearest neighbor interaction /eff estimated to be 1.1 
K.^ The ferromagnetic interaction in combination with the 
strong single- site anisotropy stabilizes the local spin arrange- 
ment in such a way that the ground state configuration has 
two spins pointing out of the tetrahedron and two spins point- 
ing into the tetrahedron, (the so-called 'two-in two-out' state 
or 'ice rules')."^ For every tetrahedron, there are six possible 
combinations of spins under the two-in two-out rule reflecting 
the global cubic symmetry. In fact, Dy2Ti207 is found to show 
the residual ground state entropy of 1.68 J/(K mole), which is 
numerically in good agreement with the Pauling's entropy for 
water ice.^'^ The spin ice state is attained by the development 
of a short-ranged spin correlations. 

Recent theoretical studies reveal that the excitations from 
the ground state can be described by the creation of magnetic 
monopoles.^ Violating the ice rules by making a spin flip on 
the ground state configuration leads to a pair of point-like de- 
fects in the tetrahedra that contain the spins. In fact, various 
experimental results are explained by the excitation of mag- 
netic monopoles.^"^^ 

The spin dynamics of Dy2Ti207 revealed by AC sus- 
ceptibility xac measurements down to 0.6 K show a quite 
unique behavior, which can be roughly described by three 
regimes. Above ~ 10 K, the temperature dependence of 
the relaxation time riT) is eff'ectively explained on the basis 
of the Arrhenius law t{T) = Tth^^viE^IT), with an energy 
barrier = 220 K.^^ In the intermediate regime, in the tem- 
perature range of 2-10 K, t{T) is almost constant, and be- 
low 2 K, t{T) increases again and reaches ~ 1 s at 0.75 K.^^ 
A diff'erence between the zero field cooled (ZFC) and field 



cooled (FC) DC magnetization is observed below 0.7 K.^^'^^ 
These dynamics of spin ice may be understood by two ori- 
gins. At high temperature above ~ 10 K, the dynamics is due 
to single ion process and mixing with excited states. 
corresponds to the energy of the excited CEF levels which 
results in the Ising-like anisotropy. At low temperature be- 
low ~ 10 K, the dynamics comes from the creation or anni- 
hilation of magnetic monopoles, and their diff'usion. Within 
the nearest neighbor spin ice model, the energy cost of a sin- 
gle spin flip is 4/eff • However, the experimental data of t{T) 
can be well fit by an Arrhenius law with an energy barrier 
of E^ = 2/eff ~ 2.2 K over the restricted temperature range 
between 2.5 and 5 K.^^'^"^ In a similar way, the low temper- 
ature increase of t(T) can also be fitted by an Arrhenius law 
but with Ep = 6/eff ~ 6.6 K in the temperature range 2 K 
down to 0.75 K. This behavior has been explained by tak- 
ing into account the long range Coulomb interaction between 
monopoles. In addition, a detailed analysis of diff'usive mo- 
tion of monopoles using classical Monte Cairo simulation can 
reproduce the t(T) dependence down to 0.75 K.^^' 

Below 0.7 K, the spin dynamics were probed by magne- 
tocaloric measurements,^^ but the large applied fields make 
the comparison with;^AC data quite difficult. Short time relax- 
ation of the magnetization has been measured down to 0.36 K 
and could be accounted for within a chemical kinetic model. 
To get a deeper insight in the spin dynamics of Dy2Ti207, we 
have performed xac measurements and DC relaxation mea- 
surements down to 0.08 K. In this letter, we show that the 
magnetic relaxation curves as well as the xac data is well 
described by using two relaxation times down to very low 
temperature. Furthermore, we show that t(T) flattens at very 
low temperature, thus deviating from the suggested Coulomb 
model. 

Single crystals of Dy2Ti207 were prepared by the floating- 
zone method, as reported previously. ^^;^ac and DC magnetic 
relaxation along the [111] direction were measured using a 
SQUID magnetometer developed at the Institut Neel, CNRS, 
Grenoble. Two samples were used to ensure the experimental 
results; the size of the samples (#1 and #2) were 0.7 x 2.2 x 2.2 
mm^ and 0.9 x 1.9 x 3.8 mm^, respectively. No meaningful 
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sample dependence is observed. In this letter, the data of ;^ac 
for sample #1 and DC magnetic relaxation measurements for 
sample #2 are shown; the data of DC magnetization at 0.08 
K is shown for sample #1. The;^AC measurements were per- 
formed using AC magnetic fields < 1.4 Oe. 

In order to minimize the demagnetizing field corrections, 
the [111] direction was oriented along the sample plane. The 
demagnetization correction for the data of Xhc is performed 
using the factor N = 0.197 x 4;: (CGS units) that was calcu- 
lated with the analytical form for a rectangular prism. The 
relaxation of the magnetization was made using the following 
protocol; (i) a weak DC magnetic field of 10 or 5 Oe is ap- 
plied, (ii) the sample is warmed up to 0.90 K for a wait period 
of 30 s. (iii) the heater power was cut, and the sample was 
subsequently rapidly cooled down, dropping below 0.4 K in 
less than 10 s and settling down to below 0. 1 K after 600 s (iv) 
the heater power was restored and regulated and stabilized at 
the target temperature for an additional 600 s. (v) Finally the 
applied field is cut (in a time ^ < 0.1 s) and the relaxation of 
the DC magnetization M(t) is measured. 

Clearly the most important parameter for the study of the 
dynamic behavior is the relaxation time t(T), which can be 
deduced from;^AC as well as the relaxation of the magnetiza- 
tion. In the case of Debye relaxation with a single dispersion, 
the imaginary part of the AC susceptibility x'' is Lorentzian, 
which is a symmetric function of log/, and a maximum oc- 
curs when the measurement time l/27if is equal to the re- 
laxation time T. In addition, the magnetization is expected to 
decay as an exponential function of time, with a characteris- 
tic time equal to r. The frequency dependence of xac in the 
temperature range from 0.5 to 1.9 K is shown in Fig. 1. Our 
data are in qualitative agreement with a previous study made 
on a poly crystalline sample between 0.8 and 1.8 K}^'^^ The 
position of the peaks in;^''(/) shift to lower frequency with 
decreasing temperature, implying that r increases on cooling 
as expected for a thermally activated process. The shape of 
the observed ;^''(/) curve is very close to a Lorentzian sym- 
metric function so that the main part of the dissipation can 
be attributed to a single relaxation time. However the curves 
do show a small but clearly discernible shoulder on their high 
frequency sides, indicating the existence of a second relax- 
ation time, or perhaps a distribution of relaxation times. Note 
that a similar feature in;^AC was observed for the high temper- 
ature data (T > 14 K). The data have been analyzed using the 
Davidson-Cole formulation with a cutoff' time, as previously 
pointed out.^^'^^ 

These observations are in agreement with the analysis of 
M(t) curves shown in Fig. 2 below 0.60 K. They were ob- 
tained (as described above) by the FC condition and then 
quenching the field at the measurement temperature. As ob- 
served in;^AC, the relaxation times increase at lower temper- 
ature. Note that, at 0.08 K, the magnetization is reduced by 
only 0.35 % after 16.7 h (See Fig. 2(b)), showing that slow 
relaxation still occurs at this temperature. 

To account for the whole set of measurements, (xac be- 
tween 0.50 and 1.9 K and M(t) below 0.60 K), we have fitted 
the data with a 2-t Debye model as follows: 

27r/TL , 27r/Ts 



X\f)=XL 



l+(27r/TL)2 ^^'l + (2;r/Ts)2' 



(1) 
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Fig. 1. (a)(Color online) Frequency dependence of x\f) of Dy 2X1267 at 
0.50, 0.55, 0.60, 0.70, 0.80 0.90. 1.0, 1.1, 1.2, 1.3, 1.5, 1.7, and 1.9 K. 
(b)(Color online) Frequency dependence of x" if) of Dy2Ti207 at 0.50, 0.55, 
0.60, 0.70, 0.80 0.90. 1.0, 1.1, 1.2, 1.3, 1.5, 1.7, and 1.9 K. 
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Fig. 2. (Color online) Relaxation of magnetization due to field quench after 
field cooling (FC) at (a) 0.25, 0.30, 0.35, 0.40, 0.45, 0.50, 0.55, and 0.60 K 
and (b) 0.20, 0.15, and 0.08 K. The applied magnetic fields in FC //pc for the 
data at 0.15-0.60 and 0.08 K are 10 and 5 Oe, respectively. 



t t 
M(t) = MLexp( ) Msexp( ) Mq, 



(2) 



where ;^L,;ts, ^s, and Mq are fitting parameters, and tl 
and Ts are the short and long relaxation times. Similar shape 
of ;^AC data in spin ice Ho2Ti207 is obtained in the tempera- 
ture range of 0.5-1.3 K?^ Figs. 3(a) and (b) show the fitting 
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results foryX/) at 0.80 K and M{t) at 0.45 K, respectively. 
The fits are rather good over a large frequency or time range 
although they are not perfect and small diff'erences can be seen 
at higher frequency and at very short or very long times. How- 
ever, the fitting results show that the bulk of the relaxation can 
be reasonably described by two relaxation times. In addition, 
the 2-T model allows to fit both AC and DC data in a wide 
temperature range (0.25-1.9 K), and thus to get an estimation 
of TL(r) and Ts(r) (See Fig. 4(a)). Below 0.25 K, the relax- 
ation of the magnetization is so slow that we cannot fit the data 
within the 2-t model. However, we can estimate the average 
relaxation time Tav using a simple exponential relaxation (See 
Fig. 4(a) inset). Note that the tl and ts at 0.55 and 0.60 K 
obtained from M{t) are in agreement with those from;^''(/). 
which shows that the 2-t model is consistent. 




Fig. 3. (Color online) Fitting results for (d.)x"{f) at 0.80 K and (b) Mit) at 
0.45 K. Solid lines (red) show the fitting curves with ixi^^Xs^^i^^ ^s) = (24.06, 
5.594, 0.1088, 0.02110) fory\f) and (Ml, Ms, Mq, tl, ts) = (119.6, 97.60, 
3.254, 380.0, 79.72) for M(t), respectively. Broken and dotted lines show the 
contribution from rs and tl, respectively. For comparison, dash-dotted lines 
(black) show the fitting curves by a single r model, which can not reproduce 
the data. 



It is interesting to speculate on the origin of two relaxation 
times in this system. To this end, we note that both relax- 
ation times have a similar temperature dependence, at least 
down to 0.3 K, but their relative contributions to x'\f) 
Mit) depend on temperature. This can be seen from the pref- 
actors;^L, X^^ and Ms in Figs. 4(b) and (c).^^ This ob- 
servation allows to conclude that the two r do not originate 
from the existence of two diff'erent sites in the pyrochlore lat- 
tice in [1 1 1] applied field.^^ A recent theoretical study on the 
basis of magnetic monopoles proposes that two types of mag- 
netic monopoles (free and bound magnetic monopoles) are 
created when a dipolar spin ice is quenched from high tem- 
perature into the spin ice regime.. In this model, the free 
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Fig. 4. (a) (Color online) Temperature dependence of tl and rs of 
Dy2Ti207. Solid and broken lines show the Arrhenius law with Ep = 9.2 
and 6.3 K, respectively. Inset shows the temperature dependence of Tav, tl 
and Ts down to 0.08 K.(b)(Color online) Temperature dependence of fitting 
parameters and;^s- (c)(Color online) Temperature dependence of Ml, Ms, 
and Mq. 



takes place when two monopoles meet elsewhere, whereas the 
bound magnetic monopoles have a longer lifetime due to a dy- 
namical arrest. However the calculations were made in zero 
field,^^ with a very fast cooling rate, presumably much faster 
than the experiment (10 s to 0.4K, several minutes to 0.1 K), 
which makes it difficult to compare directly. Notwithstanding, 
the two relaxation times in our 2-t model may correspond to 
the lifetime of these free and bound monopoles. 

Below 2 K, Tl and ts have almost the same temperature 
dependence. They increase on cooling with the largest rate 
of change between 0.5 and 1 K. Surprisingly, the increase of 
both Tl and ts tends to be suppressed below 0.5 K but the 
relaxation times continue to increase although more slowly, 
at least down to 0.08 K, reaching times as long as 200 days at 
0.08 K. 

To get an insight for the energy involved in these relaxation 
process we can describe the t(T) dependence with a Arrhe- 
nius law: t(T) = Toexp(£'p/r) in the temperature range 0.5-1 
K. From the temperature dependence of tl and ts, we get an 
energy barrier ^ 9.2 K (8 ~ 9/eff) with tql = 1.35 x 10"^ 
s and Tos = 2.3 x 10"^ s, which is larger than the previous 
result Ep ~ 6.6 K in the temperature range 0.75-2 K.^^'^^ 
This Arrhenius fit describes only a narrow temperature range. 
A description taking explicitly into account the long range 
Coulomb interaction between the monopoles is desirable to 
explain the dynamics below 0.75 K. 

A possible origin for the suppression of the rate of increase 
of t(T) below 0.5 K could be a reduction of Ep at low tem- 
perature. The reduction of the barrier comes about because 
only single defects are allowed in the low temperature limit. 
Then, for Tav below 0.2 K, assuming the Arrhenius law, a 
roush value of Er^ ~ 0.5 K is obtained. Obviouslv. this value 
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is smaller than 2/eff expected for a single defect in the near- 
est neighbor spin ice model. Therefore, it is hard to describe 
the dynamics at least below 0.2 K by the Arrhenius law. As 
another possible origin is quantum dynamic (QD) effect be- 
tween two spins. This QD effect results from the anisotropic 
superexchange coupling between neighboring Kramers dou- 
blets; the QD effect is intrinsic to pyrochlore rare-earth oxides 
and very significant in Pr pyrochlore oxides. The QD of 
Dirac string derived from the spin-flip exchange interaction 
may suppress the increasing of t(T) at very low temperature. 
This effect is expected to be smaller in the Dy2Ti207, due to 
the large spin and the strong anisotropy. A quantum mechan- 
ical treatment of the spins will be necessary to clarifying the 
dynamics due to QD at very low temperature. Note that a sim- 
ilar QD effect is observed in the nuclear quadrupole resonance 
measurement of Dy2Ti207; the increasing of spin-lattice re- 
laxation time Ti is suppressed below 0.5 K and Ti becomes 
almost constant below 0.3 K.^^ These results suggest a QD 
effect in spin ice. 

It is worth noting that our results may be complementary to 
the recent experiments that report a relaxing magnetic current 
due to monopoles.^^ These relaxation measurements were 
made after ZFC and then applying a field pulse during 5 to 
60 s. They focused on the short time relaxation (t < 180 s). 
The latter could be described using a chemical kinetic model, 
involving dissociation and recombination of the monopoles. 
The so-obtained energy barriers at 0.36 K, are 4.5 K for the 
free monopoles, and 3.3 K for the bound monopoles. On the 
contrary, our relaxation curves focus on the way the system 
reaches equilibrium at long times. As seen above, our rough 
2-T model tends to show that the long time relaxation may also 
imply free and bound monopoles. However, fitting the Arrhe- 
nius law, over a very restricted temperature range centered on 
0.36 K gives energy barriers of approximately 4 and 4.7 K 
for the short and long relaxation process respectively, and a 
very slow tq of approximately 10"^ s for both. Although the 
concentration of each type of monopoles seem to be very dif- 
ferent in our experiment, it is not surprising since in our mea- 
surements the field was applied at "high" temperature, thus 
generating different excitations than the ones created at low 
temperature. 

In summary, the xac and M(t) measurements due to field 
quench after FC along [111] direction have been performed in 
order to estimate t(T) in a spin ice compound Dy2Ti207 down 
to 0.08 K. We found that the dynamics of spin ice below 2 K 
is well described by 2-t model with ts and tl. From the view 
of magnetic monopole picture, the existence of two relaxation 
times corresponds to two types of magnetic monopoles (free 
and bound monopoles). Below 0.5 K, both ts and tl show 
the clear deviation from the thermal activated dynamics with 
£'p ~ 9 K. The increasing of the relaxation times are sup- 
pressed below 0.5 K, suggesting a QD effect. 
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